Repressive epigenetic mechanisms, such as the H3K27me3 histone modification, were predicted to affect muscle gene expression and its mineral content in Nelore cattle.

Epigenetic repression has been linked to the regulation of different cell states. In this study, we focus on the influence of this repression, mainly by H3K27me3, over gene expression in muscle cells, which may affect mineral content, a phenotype that is relevant to muscle function and beef quality. Based on the inverse relationship between H3K27me3 and gene expression (i.e., epigenetic repression) and on contrasting sample groups, we computationally predicted regulatory genes that affect muscle mineral content. To this end, we applied the TRIAGE predictive method followed by a rank product analysis. This methodology can predict regulatory genes that might be affected by repressive epigenetic regulation related to mineral concentration. Annotation of orthologous genes, between human and bovine, enabled our investigation of gene expression in the Longissimus thoracis muscle of Bos indicus cattle. The animals under study had a contrasting mineral content in their muscle cells. We identified candidate regulatory genes influenced by repressive epigenetic mechanisms, linking histone modification to mineral content in beef samples. The discovered candidate genes take part in multiple biological pathways, i.e., impulse transmission, cell signalling, immunological, and developmental pathways. Some of these genes were previously associated with mineral content or regulatory mechanisms. Our findings indicate that epigenetic repression can partially explain the gene expression profiles observed in muscle samples with contrasting mineral content through the candidate regulators here identified.

Bivariate GWAS reveals pleiotropic regions among feed efficiency and beef quality-related traits in Nelore cattle.

Feed-efficient cattle selection is among the most leading solutions to reduce cost for beef cattle production. However, technical difficulties in measuring feed efficiency traits had limited the application in livestock. Here, we performed a Bivariate Genome-Wide Association Study (Bi-GWAS) and presented candidate biological mechanisms underlying the association between feed efficiency and meat quality traits in a half-sibling design with 353 Nelore steers derived from 34 unrelated sires. A total of 13 Quantitative Trait Loci (QTL) were found explaining part of the phenotypic variations. An important transcription factor of adipogenesis in cattle, the TAL1 (rs133408775) gene located on BTA3 was associated with intramuscular fat and average daily gain (IMF-ADG), and a region located on BTA20, close to CD180 and MAST4 genes, both related to fat accumulation. We observed a low positive genetic correlation between IMF-ADG (r = 0.30 ± 0.0686), indicating that it may respond to selection in the same direction. Our findings contributed to clarifying the pleiotropic modulation of the complex traits, indicating new QTLs for bovine genetic improvement.

Prune homolog 2 with BCH domain (PRUNE2) gene expression is associated with feed efficiency-related traits in Nelore steers.

Animal feeding is a critical factor in increasing producer profitability. Improving feed efficiency can help reduce feeding costs and reduce the environmental impact of beef production. Candidate genes previously identified for this trait in differential gene expression studies (e.g., case-control studies) have not examined continuous gene-phenotype variation, which is a limitation. The aim of this study was to investigate the association between the expression of five candidate genes in the liver, measured by quantitative real-time PCR and feed-related traits. We adopted a linear mixed model to associate liver gene expression from 52 Nelore steers with the following production traits: average daily gain (ADG), body weight (BW), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), Kleiber index (KI), metabolic body weight (MBW), residual feed intake (RFI), and relative growth ratio (RGR). The total expression of the prune homolog 2 (PRUNE2) gene was significantly associated with DMI, FCR, FE, and RFI (P < 0.05). Furthermore, we have identified a new transcript of PRUNE2 (TCONS_00027692, GenBank MZ041267) that was inversely correlated with FCR and FE (P < 0.05), in contrast to the originally identified PRUNE2 transcript. The cytochrome P450 subfamily 2B (CYP2B6), early growth response protein 1 (EGR1), collagen type I alpha 1 chain (COL1A1), and connective tissue growth factor (CTGF) genes were not associated with any feed efficiency-related traits (P > 0.05). The findings reported herein suggest that PRUNE2 expression levels affects feed efficiency-related traits variation in Nelore steers.

DNA methylation may affect beef tenderness through signal transduction in Bos indicus.

BACKGROUND: Beef tenderness is a complex trait of economic importance for the beef industry. Understanding the epigenetic mechanisms underlying this trait may help improve the accuracy of breeding programs. However, little is known about epigenetic effects on Bos taurus muscle and their implications in tenderness, and no studies have been conducted in Bos indicus. RESULTS: Comparing methylation profile of Bos indicus skeletal muscle with contrasting beef tenderness at 14 days after slaughter, we identified differentially methylated cytosines and regions associated with this trait. Interestingly, muscle that became tender beef had higher levels of hypermethylation compared to the tough group. Enrichment analysis of predicted target genes suggested that differences in methylation between tender and tough beef may affect signal transduction pathways, among which G protein signaling was a key pathway. In addition, different methylation levels were found associated with expression levels of GNAS, PDE4B, EPCAM and EBF3 genes. The differentially methylated elements correlated with EBF3 and GNAS genes overlapped CpG islands and regulatory elements. GNAS, a complex imprinted gene, has a key role on G protein signaling pathways. Moreover, both G protein signaling pathway and the EBF3 gene regulate muscle homeostasis, relaxation, and muscle cell-specificity. CONCLUSIONS: We present differentially methylated loci that may be of interest to decipher the epigenetic mechanisms affecting tenderness. Supported by the previous knowledge about regulatory elements and gene function, the methylation data suggests EBF3 and GNAS as potential candidate genes and G protein signaling as potential candidate pathway associated with beef tenderness via methylation.

Muscle allele-specific expression QTLs may affect meat quality traits in Bos indicus.

Single nucleotide polymorphisms (SNPs) located in transcript sequences showing allele-specific expression (ASE SNPs) were previously identified in the Longissimus thoracis muscle of a Nelore (Bos indicus) population consisting of 190 steers. Given that the allele-specific expression pattern may result from cis-regulatory SNPs, called allele-specific expression quantitative trait loci (aseQTLs), in this study, we searched for aseQTLs in a window of 1 Mb upstream and downstream from each ASE SNP. After this initial analysis, aiming to investigate variants with a potential regulatory role, we further screened our aseQTL data for sequence similarity with transcription factor binding sites and microRNA (miRNA) binding sites. These aseQTLs were overlapped with methylation data from reduced representation bisulfite sequencing (RRBS) obtained from 12 animals of the same population. We identified 1134 aseQTLs associated with 126 different ASE SNPs. For 215 aseQTLs, one allele potentially affected the affinity of a muscle-expressed transcription factor to its binding site. 162 aseQTLs were predicted to affect 149 miRNA binding sites, from which 114 miRNAs were expressed in muscle. Also, 16 aseQTLs were methylated in our population. Integration of aseQTL with GWAS data revealed enrichment for traits such as meat tenderness, ribeye area, and intramuscular fat . To our knowledge, this is the first report of aseQTLs identification in bovine muscle. Our findings indicate that various cis-regulatory and epigenetic mechanisms can affect multiple variants to modulate the allelic expression. Some of the potential regulatory variants described here were associated with the expression pattern of genes related to interesting phenotypes for livestock. Thus, these variants might be useful for the comprehension of the genetic control of these phenotypes.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes.

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

Muscle transcriptome analysis reveals genes and metabolic pathways related to mineral concentration in Bos indicus.

Mineral content affects the biological processes underlying beef quality. Muscle mineral concentration depends not only on intake-outtake balance and muscle type, but also on age, environment, breed, and genetic factors. To unveil the genetic factors involved in muscle mineral concentration, we applied a pairwise differential gene expression analysis in groups of Nelore steers genetically divergent for nine different mineral concentrations. Here, based on significant expression differences between contrasting groups, we presented candidate genes for the genetic regulation of mineral concentration in muscle. Functional enrichment and protein-protein interaction network analyses were carried out to search for gene regulatory processes concerning each mineral. The core genetic regulation for all minerals studied, except Zn, seems to rest on interactions between components of the extracellular matrix. Regulation of adipogenesis-related pathways was also significant in our results. Antagonistic patterns of gene expression for fatty acid metabolism-related genes may explain the Cu and Zn antagonistic effect on fatty acid accumulation. Our results shed light on the role of these minerals on cell function.

Oxidative stress and antioxidant status in beta-thalassemia heterozygotes.

BACKGROUND: Several studies have evaluated the oxidant and antioxidant status of thalassemia patients but most focused mainly on the severe and intermediate states of the disease. Moreover, the oxidative status has not been evaluated for the different beta-thalassemia mutations. OBJECTIVE: To evaluate lipid peroxidation and Trolox equivalent antioxidant capacity in relation to serum iron and ferritin in beta thalassemia resulting from two different mutations (CD39 and IVS-I-110) compared to individuals without beta-thalassemia. METHODS: One hundred and thirty subjects were studied, including 49 who were heterozygous for beta-thalassemia and 81 controls. Blood samples were subjected to screening tests for hemoglobin. Allele-specific polymerase chain reaction was used to confirm mutations for beta-thalassemia, an analysis of thiobarbituric acid reactive species was used to determine lipid peroxidation, and Trolox equivalent antioxidant capacity evaluations were performed. The heterozygous beta-thalassemia group was also evaluated for serum iron and ferritin status. RESULTS: Thiobarbituric acid reactive species (486.24 ± 119.64 ng/mL) and Trolox equivalent antioxidant capacity values (2.23 ± 0.11 mM/L) were higher in beta-thalassemia heterozygotes compared to controls (260.86 ± 92.40 ng/mL and 2.12 ± 0.10 mM/L, respectively; p-value < 0.01). Increased thiobarbituric acid reactive species values were observed in subjects with the CD39 mutation compared with those with the IVS-I-110 mutation (529.94 ± 115.60 ng/mL and 453.39 ± 121.10 ng/mL, respectively; p-value = 0.04). However, average Trolox equivalent antioxidant capacity values were similar for both mutations (2.20 ± 0.08 mM/L and 2.23 ± 0.12 mM/L, respectively; p-value = 0.39). There was no influence of serum iron and ferritin levels on thiobarbituric acid reactive species and Trolox equivalent antioxidant capacity values. CONCLUSION: This study shows an increase of oxidative stress and antioxidant capacity in beta-thalassemia heterozygotes, mainly in carriers of the CD39 mutation.

Oxidative stress and antioxidant status in beta-thalassemia heterozygotes

Background: Several studies have evaluated the oxidant and antioxidant status of thalassemia patients but most focused mainly on the severe and intermediate states of the disease. Moreover, the oxidative status has not been evaluated for the different beta-thalassemia mutations. Objective: To evaluate lipid peroxidation and Trolox equivalent antioxidant capacity in relation to serum iron and ferritin in beta thalassemia resulting from two different mutations (CD39 and IVS-I-110) compared to individuals without beta-thalassemia. Methods: One hundred and thirty subjects were studied, including 49 who were heterozygous for beta-thalassemia and 81 controls. Blood samples were subjected to screening tests for hemoglobin. Allele-specific polymerase chain reaction was used to confirm mutations for beta-thalassemia, an analysis of thiobarbituric acid reactive species was used to determine lipid peroxidation, and Trolox equivalent antioxidant capacity evaluations were performed. The heterozygous beta-thalassemia group was also evaluated for serum iron and ferritin status. Results: Thiobarbituric acid reactive species (486.24 ± 119.64 ng/mL) and Trolox equivalent antioxidant capacity values (2.23 ± 0.11 mM/L) were higher in beta-thalassemia heterozygotes compared to controls (260.86 ± 92.40 ng/mL and 2.12 ± 0.10 mM/L, respectively; p-value < 0.01). Increased thiobarbituric acid reactive species values were observed in subjects with the CD39 mutation compared with those with the IVS-I-110 mutation (529.94 ± 115.60 ng/mL and 453.39 ± 121.10 ng/mL, respectively; p-value = 0.04). However, average Trolox equivalent antioxidant capacity values were similar for both mutations (2.20 ± 0.08 mM/L and 2.23 ± 0.12 mM/L, respectively; p-value = 0.39). There was no influence of serum iron and ferritin levels on thiobarbituric acid reactive species and Trolox equivalent antioxidant capacity values. Conclusion: ...

F200Y polymorphism in the ß-tubulin gene in field isolates of Haemonchus contortus and risk factors of sheep flock management practices related to anthelmintic resistance.

Haemonchus contortus is the most prevalent and pathogenic nematode of sheep in tropical areas. The objectives of this study were to assess the frequency of the F200Y polymorphism on the ß-tubulin gene in third-stage larvae of H. contortus from 33 sheep flocks in São Paulo state, Brazil, and to associate this frequency to risk factors based on farm management practices. The resistance allele frequency varied from 9 to 74%, and the resistance genotype frequency varied from 0 to 66.7%. Resistance genotype frequencies higher than 40% were associated with multiple risk factors - new sheep farming enterprises, the absence of farm records, the use of Dorper and Suffolk breeds, rotational grazing, the lack of wetlands on farms, pasture sharing with cattle or horses, frequent incorporation of animals into the flock, semi-intensive farming systems, whole-flock treatment, failure to use the FAMACHA method, lack of the dose-and-move practice, anthelmintic rotation after each application, visual estimation of animal weight for treatment, and lack of drug combination use. It can be concluded that genotyping the F200Y polymorphism can be used to monitor the resistance in sheep flocks and the knowledge of management strategies at the farm level is important to identify drug resistance related factors.

Polymorphism and parent-of-origin effects on gene expression of CAST, leptin and DGAT1 in cattle.

This study aimed to investigate differential allele expression (DAE) and polymorphism and parent-of-origin effects on expression of genes related to beef traits. CAST, related to meat tenderness, and DGAT1 and leptin, related to fat deposition, were evaluated. In bovine fetal tissues CAST was expressed twice as much (P < 0.05) in muscle of homozygous GG than in heterozygous AG. Leptin was expressed about one-tenth as much (P < 0.05) in heterozygous TpCm (allele T of paternal origin and allele C of maternal origin) than in homozygous CC. No DAE was observed. The evidence of polymorphism effect on expression of CAST and parent-of-origin effect on leptin contributes to a better understanding of events controlling the expression of genes of economic interest in cattle. Furthermore, if the parent-of-origin effects observed in fetal tissues are confirmed in adult tissues and associated to phenotypic variation, this parental origin criterion may be considered in marker-assisted selection of beef traits.